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Abstract 



We consider the open XXZ quantum spin chain with nondiagonal boundary terms. For bulk 
anisotropy value rj = p = 1 , 2 , . . ., we propose an exact {p + l)-order functional relation 
for the transfer matrix, which imphes Bethe-Ansatz-hke equations for the corresponding 
eigenvalues. The key observation is that the fused spin-^^^ transfer matrix can be expressed 
in terms of a lower-spin transfer matrix, resulting in the truncation of the fusion hierarchy. 



1 Introduction 



A long outstanding problem has been to solve the open XXZ quantum spin chain with 
nondiagonal boundary terms, defined by the Hamiltonian f!], |, |^ 



^ N-l 



n=l 

+ sinhr/(^coth^_ai^ + — ^a^-cothe+a^-— ^a^jj, (1.1) 

where cr^ , o"^ , cr^ are the standard Pauli matrices, ri is the bulk anisotropy parameter, ^± , k± 
are arbitrary boundary parameters, and is the number of spins. Solving this problem (e.g., 
determining the Bethe Ansatz equations) is a crucial step in formulating the thermodynamics 
of the spin chain and of the boundary sine-Gordon model. Moreover, this problem has 
important applications in condensed matter physics and statistical mechanics. 

A fundamental difficulty is that, in contrast to the special case of diagonal boundary 
terms (i.e., k± = 0) considered in 0, |l|, a simple pseudovacuum state does not exist. Hence, 
most of the techniques which have been developed to solve integrable models cannot be 
applied. Moreover, it is not yet clear how to implement the few techniques (such as Baxter's 
T — Q approach |P or the generalized algebraic Bethe Ansatz p, |^) which do not rely on a 
pseudovacuum state. 

We report here some progress on this problem. Namely, for bulk anisotropy value 

ZTT 

77=—-, p=l,2,..., (1.2) 
p+l 

(and hence q = e'' is a root of unity, satisfying q^'^^ = —1), we propose an exact (p+ l)-order 
functional relation for the fundamental transfer matrix, which implies Bethe-Ansatz-like 
equations for the corresponding eigenvalues. The key observation is that the fused transfer 
matrices t^^\u), which are constructed with a spin-j auxiliary space, satisfy the identity Q 



t(V)(„) = a(M) t^'^\u + 7]) + /3{u)I , (1.3) 

where t^^^ = I (identity matrix), and a{u) , (3{u) are scalar functions. That is, the spin-^^i 

transfer matrix can be expressed in terms of a lower-spin transfer matrix, resulting in the 

truncation of the fusion hierarchy. We have verified this result explicitly for p = 1,2,3, 

and we conjecture that it is true for all positive integer values of p. The simplest case 

^This is distinct from the observation due to Belavin et al. |0, ||] that, for the special case of quantum- 
group symmetry (i.e., k± = 0, ^± ^ cxd), the fused transfer matrix t'^i\u) vanishes after quantum group 
reduction. 
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p = 1, which corresponds to the XX chain, has recently been analyzed in 0. Similar higher- 
order functional relations have been obtained for the closed (periodic boundary conditions) 



8-vertex model by Baxter |]I0[ using a different method. 

The outline of this article is as follows. In Section ^ we review the construction of the 
fundamental {j = |) transfer matrix which contains the Hamiltonian (|1 . 1| ) . In Section ^ 
we briefly review the so-called fusion procedure [^-[|T^ and the construction of the higher- 
spin transfer matrices, which obey an infinite fusion hierarchy. In Section ^, we obtain the 
relation ( p..3|) , to which we refer as the "truncation identity," since it serves to truncate the 
fusion hierarchy. In Section ^ we present the exact functional relations which are obeyed 
by the fundamental transfer matrix. The corresponding sets of Bethe Ansatz equations are 
given in Section We conclude with a brief discussion of our results in Section ^. 



2 Fundamental transfer matrix 

We recall [|1| that the transfer matrix for an open chain is made from two basic building 
blocks, called R (bulk) and K (boundary) matrices. 

An R matrix is a solution of the Yang-Baxter equation 



Ruiu - v) R^:^(u) Ry.^iv) = R-?.:i(v) R^:^(u) R^'?iu - v) . 



1-131 



•-23 



231 



•-131 



1-121 



(2.1) 



(See, e.g., [13, |T^, [T9[-) For the XXZ spin chain, the R matrix is the 4x4 matrix 

\ 

D sinh ?/ sinh n fl 

R(u) 



( sinh(u + 77) 

sinhu smhrj 

sinh?7 sinh-u 

smh{u + 7]) J 



(2.2) 



\ 

where rj is the anisotropy parameter. This R matrix has the symmetry properties 

Ri2{u) = VuRMVu = RuiuY^'' , 



(2.3) 



where V12 is the permutation matrix and t denotes transpose. Moreover, it satisfies the 
unitarity relation 



Ri2{u) Ri2{—u) = C{u)I , ({u) = — sinh(M -|- rj) sinh(tt — rj) , 



(2.4) 



and the crossing relation 



Ru{u) = V,R,2{~u - riY^V, , V = -tcry . 



(2.5) 
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Finally, it has the periodicity property 

Ruiu + in) = -alRi2{u)al = -alRi2iu)al . (2.6) 

The matrix K~ (u) is a solution of the boundary Yang-Baxter equation 



Ru{u - v) K-{u) R2i{u + v) K~{v) = K^{v) Ru{u + v) K^{u) R2i{u - v) . (2.7) 
We consider here the following 2x2 matrix 0, ^ 

^_^^^^/sinh(e- + ..) /.-sinh2. \ ^^.8) 
\ K_sinh2u sinh((^_ — m) I 

which evidently depends on two boundary parameters , We set the matrix K~^{u) to 
be K~{—u — r]) with ((^_ , /t_) replaced by ((^+ , /€+); i.e., 

^ f - smh{u + sinh(2u + 2r/) \ 

I — K+ sinh(2u + 27]) sinh(u + r] + ^^) J 

The K matrices have the periodicity property 

K^iu + in) = -a'K^{u)a\ (2.10) 

The fundamental transfer matrix t{u) for an open chain of spins is given by [jl| 

t{u) = tio K+iu) To{u) K^{u) fo{u) , (2.11) 

where tro denotes trace over the "auxiliary space" 0, and Tq{u), To (A) are so-called mon- 
odromy matrices [| 

To{u) = Ron{u) ■ ■ ■ Roi{u) , fo{u) = Rio{u)---RNo{u). (2.12) 

Indeed, Sklyanin has shown that t{u) constitutes a one-parameter commutative family of 
matrices 

[t{u) ,t{v)] = 0. (2.13) 
The Hamiltonian ( |1 . 1|) is related to the first derivative of the transfer matrix 



n = ^^^L^^ sinh^^ + iVcosh^^^ ^2.14) 

4 sinh ^_ sinh ^+ sinh r] cosh rj 2 cosh r] 



2 



As is customary, we usually suppress the "quantum-space" subscripts 1 , . . . , N. 
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The corresponding energy eigenvalues E are therefore given by 

A'(0) sinh^ ?7 + cosh^ 7] 



2 cosh 77 



(2.15) 



£j 

4 sinh ^_ sinh ^+ sinh^^~^ 77 cosh rj 

where A{u) are eigenvalues of the transfer matrix. 
The transfer matrix has the periodicity property 

t{u + in) = t{u) , (2.16) 

as follows from (|2.6|) , (|2.1CI|) . Moreover, the transfer matrix has crossing symmetry 

t{-u - 7]) = t{u) , (2.17) 



which can be proved using a generalization of the methods developed in for the spe- 
cial case of quantum-group symmetry. Finally, we note that the transfer matrix has the 
asymptotic behavior (for k± 7^ 0) 



gn(2Af+4)+r?(iV+2) 
t{u) ~ -K.K+ 1 + . . . 



for 



u — > 00 . 



(2.18) 



3 Fusion 

Our main tool is the so-called fusion technique, by which higher-dimensional representations 
can be obtained from lower-dimensional ones. The fusion technique was first developed in 
111 , 0, |13| for R matrices, and then later generalized in - for K matrices. Following 



T^ , we introduce the (undeformed) projectors 



ml 



(3.1) 



where the sum is over all permutations a = (cri , . . . dm) of (1 , . . . , m), and Va is the permu- 
tation operator in the space (8>^iC^. For instance. 



P. 



12 



123 



(I + P12), 

(I + P23^12 + P12P23 + Vu + V23 + P13) 



The fused spin-(j, ^) R matrix (j = ^ 1 



3 
' 2 ' 



is given by |T^, |T^ 



R{1...2j)2j+l{u) = 2i^l,2j+l(«)^2,2i+l(M + V) ■ ■ ■ R2j,2j+l{u + {2j - l)r])P^ 



2i 



(3.2) 



(3.3) 
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We also note 

R2j+l{1...2j)iu) = 2jR2j+l,2jiu - {2j - l)r]) . . . R2j+l,2iu - r])R2j+lA{u)P^ 2j 

= i?(i...2,)2,+i(«-(2j-l)r/). (3.4) 

The fused spin-j K~ matrix is given by |T5| , |16 

^(~i...2j)(«) = PL2jK2j{u)R2j,2j-i{2u + rj)K2j_i{u + r]) 

X i?2j,2j-2(2M + 2r/)i?2i-i,2i-2(2ti + 3r))K-_2{u + 2r]) 



X ... R2j,i{2u + {2j - l)r])R2,-i,i{2u + 2j7]) . . . R2,i{2u + {Aj - 3)r]) 
X Kr(«+(2j-l)r7)F+,2,. (3.5) 
The fused spin-j matrix is given by K'^^ ^.^{—u — 2jvi) with (^_ , /?_) replaced by (^+ , «;+), 



K 



(t..2,>M =^(i...2i>(-^-2jr;) 



(3.6) 



The fused (boundary) matrices satisfy generalized (boundary) Yang-Baxter equations. 
The fused transfer matrix t^^^ [u) constructed with a spin-j auxiliary space is given by 

t^'\u) = tri...2, ir+ ..2,)(«)^(i...2.)(«)^(1..2,)(^)^(i...2.>(« + (2j - m , (3.7) 

where 



T{1...2j){u) 

r(i...2,)(w + (2j-i)77) 



R{1...2j)N{u) . . . R{1...2j)l{u) , 
R{1...2j)l{u) . ..R{l...2j)N{u) . 



The transfer matrix ( |2.11| ) corresponds to the fundamental case j = |; that is, t*^2)(-u) = t{u). 
The fused transfer matrices constitute commutative families 

[t^^\u) ,t^^\v)] = . (3.9) 

These transfer matrices also satisfy a so-called fusion hierarchy [jl5|, |1^] 

t^'\u) = C2j~i{2u + {2j - 1)7]) y-^\u) t^'^\u + {2j - 1)77) 

- ^^uH2j- 2),)Cj-2(2u ^ (2j-2),) , 

a2u + 2{23-l)r]) ^ ' ^ ' 

with = I, and j = 1 , | , . . .. The quantity A('u), the so-called quantum determinant 
El, H, is given by 



A(m) = A {K+{u)] A {K-{u)]5{T{u)}5\f{u) 



(3.11) 



where 

6{T{u)} = tTu{Pi2T,{u)T2{u + r])} = au + rj)^, 

A {K-iu)} = tri2 {Pu K^{u) Ri2{2u + r,) K^{u + ri) ] 

= — sinh 2u [smh('U + 77 + ^S) sinh('U + 77 — ^J) + sinh^(2ii + 2?])] , 

A = tri2 {P{2 Kt{u + 77) i?i2(-2u - 377) 7^1+ (^x)} 



A{K-(-n-277)} (3.12) 

(5- 



Moreover, 



0M = nC(« + ^?7), CoM = l. (3.13) 

k=l 

4 Truncation identity 

We now proceed to formulate the important identity (|1.3|) , which serves to truncate the 
fusion hierarchy (|3.1CI|) . To this end, we first derive separate "truncation" identities for the 
R and K matrices. 

4.1 R matrix truncation 

We recall that, in addition to the fusion approach described above, there is an alternative 
construction |2^ of higher-spin R matrices based on quantum groups. Following the notation 
of [|], the spin- ( I , j) -R matrix is given by 

(sinh ( It + (^ + i/)r? ) sinh 77 F , 

^ ^ ^ ( \ \ . (4.1) 

sinh 77 E sinh \ u\{\- H)ri 

where the matrices E and F have matrix elements 

{H)mn = ij + l-n)Sm,n, m , 71 = 1 , 2 , . . . , 2 j + 1 , 



(-E)mn = ^m^m,n-l , {F)mn = (^n^m-l,n , = ^ Mq [2j + 1 - n]q , (4.2) 

and 

N, = ^f^, q = e^. (4.3) 



These matrices form a {2j + 1) -dimensional representation of the Uq{su{2)) algebra 

[H,E] = E, [H,F] = -F, [E,F] = [2H]g. (4.4) 
The corresponding spin-(j A) R matrix i?^^ j_ is then given by 

{RII./u)) = .^{u)) p,^. . (4.5) 



0/3 



We refer to these R matrices as "quantum group" {qg) R matrices in order to distinguish 
them from the fused R matrices constructed previously ( p.3| ). The two sets of R matrices 
are related as follows 



Bl...2jAi„,2j R(1...2j)2j+l{u) Al].2jBi],2j 



■2i-i 

JJ^ sinh(u + krf) 

.k=l 



Rll,^{u + {2j-l)j), (4.6) 



where ^i...2j is the matrix of (unnormalized) Clebsch-Gordon coefficients in the decomposi- 
tion of the tensor product of 2j spin-i representations into a direct sum of su{2) irreducible 
representations. Moreover, -Bi...2j is a w-independent diagonal matrix which renders sym- 



metric the matrix on the LHS of ( [4.6|) . The A and B matrices for the cases j = 1 , | , 2 are 
given in Appendix ^ 

The quantum-group R matrices have a particularly simple "truncation" property. Indeed, 
we find that for r] = the spin-(2±i ,^)R matrix takes the block-diagonal form 



^2 '2> 



^ i cosh(-u + 2)cr^ 










R 



19 

2 ' 2 -' 





u 








cosh('U + ^)cr^ y 



(4.7) 



In view of the relation ([4.6|) , we see that the corresponding fused spin-(2±i , i) R matrix 
satisfies 



2 ' 2' 



-Bi...p+iAi...p+i R{i...p+i)p+2{u) 

/ z/(M)cr^ 

fii...p-iv4i...p_i R{i...p-i)p{u + 7]) ^r^.p-i-^L.Vi 
\ 



\ 



,(4. 



where 



H{u) = C(n) , 



(4.9) 



■^It is understood that the "null" rows and columns (i.e., those with only zero matrix elements) are to be 
pruned from the LHS. We have verified this relation explicitly for j = 1 , | , 2, and we conjecture that it is 
true for all j . 



7 



and 



1 ^ 

— I I sinhf-u + krf) 



^ - ) sinh((p + l)u) . 



k=0 



(4.10) 



In obtaining this result, we have used the identity (see 1.392 in 



JJ^ sinh(M + krf) 



k=0 



p+i 



sinh((p + l)u) 



(4.11) 



As will be explained in Section below, it is actually more useful to consider the 
similarity transformation with a matrix C (instead of B), which results in a triangular 
(instead of a symmetric, block-diagonal) matrix. For later convenience, we present this 
result here: 



Cl...p+lAl...p+l /?(!.. .p+l)p+2('u) ^l...p+lC'l...p+l 



v[u]a 



\ 



Sl,..p_lAl.„p_l i?^l,..p_l)p(M + ?7) * 

\ -iy{u)a^ J 



(4.12) 



The monodromy matrices therefore obey analogous relations 

Cl...p+l^l...p+l T(^l...p+l){u) ^L.^.p+l^L.^p+l 



/i(tt 



ij{u)^F \ 

{-iy{u))^F j 



/ v(uYF 



\ 



Si.„p_iAi,..p_i r(i,..p_i)(w + {p- 1)7]) 



V 



-u{u))^F ) 



(4.13) 



where F = Y[k=i ^k- The C matrices for the cases j = 1 , | , 2 are also given in Appendix ^ 



4.2 K matrix truncation 



An explicit construction of nondiagonal higher-spin K matrices analogous to ( [4.1|) is unfor- 
tunately not yet known. Nevertheless, considerable insight can be gained by first considering 
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the diagonal case (i.e., k± = 0). The spin-j diagonal K matrix is given explicitly by (see 

(4.14) 



IJ] for the case j = 1) 

-q 

U) 



K-;{\u) = drag ( ^^(n) , ^g}(n) , . . . , k^"^ ^ ) , 



where 



2j-m 



m-2 



k^^\u) = Y\ sinh(^_ + u + It]) JJ^ sinh(^_ — u — Irj) . 

1=0 1=0 

This K matrix is related to the fused K matrix (|3.5| ) by 



(4.15) 



■2j-l I 

Yl JJsinh(2M+ (/ + k)ri) 

. 1=1 k=l 



(4.16) 



We find that the diagonal K^^^^{u) with j = has the following "truncation" property 



for rj 



P+i' 



K (u) = sinh(,^_ + u) sinh(,^_ — u) 

^ 2 > 



\2) si 



X 



sinh((p+l)(g-+^)) 



sinh(5_ +u) sinh(5_ —u) 














\ 



\ 



• (4.17) 



n n / i\v sinh((p+i)(g--»)) 

^ ^ \ 2> smh(5_+«)sinh(5_--u) / 



It follows that for rj = the diagonal fused spin-^ii K matrix satisfies the truncation 
identity 

-1 



u_{u) \ 









where 



2p 



(4.18) 



A{K-(u-r])} TT , , 

^ ^ n sinh 2m + kri) 

sinh2M-2?7 11 ^ 



A{K-(n - T])} sinh^ (2(p + l)u) 



sinh 2t/ sinh(2ti + 77) sinh(2'u, — 77) sinh(2ti — 2//) 



(4.19) 



^Footnote || applies here too. 
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the quantum determinant A{K (u)} is given in Eq. ( 3.12 ) 

p I 



1 



P-[U) 



Y[Ylsmh{2u + {I + k)r]) 



,1=1 k=l 



sinh((p+ 1)(^_ + «)) 



1 g|i7rp(p+2) 

— - cosh'2l((p + l)-u) sinh' 2 ^((p + l)^) 



X sinh((^»+ 1)(^_ + u)) , 



and 



z/_ Ui 



with 



(4.20) 



(4.21) 



The standard notation [x] denotes integer part of x. 

We turn now to the more general nondiagonal case {k± 7^ 0). The boundary matrices 

hke their bulk counterpart ( [4.6| ), are symmetric. However, for 77 = and j = 2±i, these 
matrices are not block-diagonaL In order to obtain a truncation identity for the full transfer 
matrix, it is desirable to have at least block-triangular matrices. 

To this end, we consider a new similarity transformation, replacing the matrix i? by a new 
matrix C which is also diagonal and w-independent. (See Appendix ^for explicit expressions 
of the C matrices for j = 1 , | , 2.) Indeed, we propose the following generalization of ( [4.18|) 
for the nondiagonal case: Q 

Ci...p+iAi...p+i A'^''^ p^-^^('u) ^i.^.p+iCi..^p+i 







,(4.22) 



where ^i-{u) and v'_{u) are given by the same expressions (|4.19| ), ( |4.21| ); but the expression 
( |4.20| ) for V-iu) is now replaced by 



cosh^^\{p + l)u) smh^^\{p + l)u) n(w;^_,/«_), (4.23) 



^This equation is formal, since CT:p — > 00 (and 0) in the limit ij as can be seen from Eqs. 



(4.26) and (4.28). However, substituting this result into the expression for the fused transfer matrix gives a 

limit. 



final result (4.30) which is finite in the j] 



P+i 
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where the function n{u ; , k) is defined by 

i e , «;) = sinh {{p + + u)) + Cp,i n"^ sinh ((p +l)u+{p+l- 2l)i) . (4.24) 



n[u : 



1=1 



We have exphcitly computed the coefficients Cp^i which appear in this function for values of 
p up to p = 5, and we find that they are consistent with the following formulas: 



Cp,l = 

1 , 

Cp,2 = 



(4.25) 



and also 05^3 = 2. It remains a challenge to determine the coefficients Cp^i for all values oi p 
and /. Moreover, 



CT-iU] 



//_(m) 



P-W 



(4.26) 



where 



p+1 



Yl Yl sinh(2M +{l + k- l)r]) 



,1=1 k=l 



P 



Y[ sinh(2?/ + 2kr]) 



. k=0 

p+3 i 



— — K^^^ cosh[Vl((p _^ sinh[Vl((p _^ 1)^) 



(4.27) 



and 



(4.28) 



Note that a 00 for r? — > ^f-. 

' p+i 

The matrices are given, in view of Eq. ( p.6| ) with 77 = and j = by 

/<-;..,«,(..) = A<L,«,(-«--)|,_,._,_,^^,.^,. (4.29) 

Hence, the "plus" quantities etc. can be readily obtained from the corresponding 

"minus" quantities etc. by making the replacements u —* —u — in and (^_ ,/€_)—*> 



4.3 Transfer matrix truncation 



We are finally in position to formulate the truncation identity for the fused transfer matrix 
t^^\u) defined in ( p.7| ). Recalling the results ( [4.13| ) for the monodromy matrices and ( [4.22| ) 
for the K matrices, we obtain (for rj = and j 



P+l ^ 
2 ' 



a[u] 



U + T]) + (3{u)l 



(4.30) 
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where 



a{u) = ii{uf^ ii^{u)iJ.^{u) , 

(3{u) = p{uf''\v_{u)p+{u) + u'_{uy^{u) 



(4.31) 



Note that the factors of a from (y^{u) and p^{u) cancel in the expression for (3{u)] and hence, 
the result is finite for r] — > 



5 Functional relations 



Combining the fusion hierarchy (p.lO|) and the truncation identity ( [4.30|) , it is straightforward 
to obtain - for any positive integer value of p - a {p + l)-order functional relation for 
the fundamental transfer matrix. We propose the following general form of the functional 
relations: 

fo{u)t{u)t{u + r]) ...t{u + pifj 

- fi {u)t{u + r])t{u + 2ri)...t{u + {p-l)ri) 

- fi{u + vi)t{u + 2r])t{u + 377) . . . t{u + prj) 

- fi{u + 2r])t{u)t{u + 3r])t{u + Ar]) . . .t{u + prj) 

- fi{u + ?)rj)t{u)t{u + vi)t{u + Arj) . . . t{u + prj) — . . . 

- fi{u + pr])t{u)t{u + rj) . . .t{u + {p - 2)rj) 
Mu)Mu + 2rj) 



+ 



+ 



+ 



+ 



+ 



+ 



foiu) 

f,(u)Mu + 3rj) 
fo{u) 

f,(u)Mu + {p-l)v) 



foiu) 

fiiu + r])fi{u + 3r]) 

/oH 
fi{u + r])fi{u + Ar]) 

foiu) 

fiju + v)fiiu + pr]) 
foiu) 



t{u + 3r])t{u + 4?7) . . . t{u + {p - I)?]) 
t{u + r])t{u + 4r])t{u + 5r]) . . . t{u + {p - l)r]) + 
t{u + r])t{u + 2r]) ... t(u + {p - 3)r]) 



t{u + 'irj)t{u + 5//) . . . t{u + prj) 
t{u + 2rj)t{u + 5ri)t{u + 677) . . . t{u + prj) + 
t{u + 2ri)t{u + 3r])... t{u + {p - 2)r]) 
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-t[u)t[u + brj)t[u + Qrj) . . . t[u + pi]) 



+ -"^ y 'U {u)t{u + 3r])t{u + 67])t{u + 7r])...t{u + pr]) + ... 





Jo{u) 






+ 2ri)h{u 


+ 577) 




fo{u) 






+ 2T,)h{u 


+ pr]) 




fo{u) 











t{u)t{u + 'ir])t{u + A7])... t{u + {p- 2)77) 

+ ... = h{u). (5.1) 

In particular, the first three functional relations are given by 

p=l: fo{u)t{u)t{u + 7])- fi{u) - fi{u + 7]) = f-i{u) , 

p = 2: foiu)t{u)t{u + r])t{u + 2ri) - fiiu)t{u + r]) - h{u + ri)t{u + 2r/) 

-h{u + 2r^)t{u) = h{u), 
p = 3: fo{u)t{u)t{u + r])t{u + 27])t{u + 3r/) - fi{u)t{u + r])t{u + 2r]) 

-fi{u + r])t{u + 2r])t{u + 3r]) - fi{u + 2r])t{u)t{u + 3r]) 
-Mu + 37])t{u)t{u + v) = h{u) . (5.2) 

Remarkably, only three distinct functions /o, /i, /s appear in the functional relations. 
The function fo{u) is given by 

p I 

1=1 k=l 

= 2Pip~i) coshP((p+ l)7/)sinhP \{p+l)u)cosh{{p+l){u- —)). (5.3) 
This function has the periodicity property 

foiu) = fo{u + v) (5.4) 
for 77 = The function fi{u) is given by 

fo{u)A{u - 77) 

/iM = ^^^^;^ ' ^^-^^ 

where C(w) is given by ( |2.4|) , and the quantum determinant A('u) is given by (|3.11| ). Finally, 
the function f3{u) is given by 

/3(n) = a{u)P{u) 
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= ^2Ztm ' cosh^^§]((p+ l)..)sinh^^+^t^l((p+ 1)^) 

+2{-lf{-ti^K+f+^ sinh2(2(p + , (5.6) 

where the function n{u;^,K,) is defined in Eq. ( [4.24|) . As aheady mentioned, we have 
computed the coefficients Cp^i which appear in this function only up to p = 5. 



We remark that, in obtaining the result (|5.1| ), we have made use of the relation 



p-2 I 

Mu) = a{u) l[l[a2u+{l + k + 2)t,) (5.7) 

1=1 k=l 



which is satisfied by the function fi{u) defined in (|5.5| ). 



6 Eigenvalues and Bethe Ansatz equations 

The functional relations which we have obtained can be used to determine the eigenvalues 
of the transfer matrix. The commutativity relation ( ^.13[ ) implies that the transfer matrix 
has eigenstates |A) which are independent of u, 

t(«)|A) = A(w)|A), (6.1) 

where A (it) are the corresponding eigenvalues. Acting on |A) with the functional relation 
( ^.1| ), we obtain the corresponding relation for the eigenvalues 

/o(u)A(u)A(u + ?7) . . . A(u + pri) - fi{u)A{u + r])A{u + 2r])... A{u + [p - l)r]) 

+ ... = f,(u). (6.2) 



Similarly, it follows from (|2.16|) and (|2.17|) that the eigenvalues have the periodicity and 
crossing properties 

A(m + iyr) = A(m) , A{-u - 7]) = A{u) . (6.3) 
Finally, ( p.l8|) implies the asymptotic behavior (for Kj- ^ 0) 

A{u) ~ 22N+1 ^ • • • M — > oo . (6.4) 

We shall assume that the eigenvalues have the form 

N 

A{u) = p Y\ sinh(M — Uj) sinh('U + rj + Uj) , (6.5) 

i=-i 
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where Uj and p are (m- independent) parameters which are to be determined. Indeed, this 
expression satisfies the periodicity and crossing properties ( |6.3|) , and it has the correct asymp- 
totic behavior ( |6.4| ) provided that we set 

p = — 8/t_K+ . (6.6) 

Evaluating the functional relation ( |0|) at the root u = uj, we obtain a set of Bethe- 
Ansatz-like equations 

-fi{uj)A{uj + r])A{uj + 27])... A{uj + {p - l)r]) 
+ /iK)/i(^^ + 2^) a(^^. + 3v)A{u, + 477) • • • A(% + (P ' ^)v) 

+ IiMM^^IIMa^Uj + n)Aiu, + Ari)Aiu, + 5ri) . . . A{u, + (p - I)?]) + . . . 
^ h(u>)h{u, + (p-l)r,) ^^^^ ^ ^^^^^^^ _^ ^ ^ ^ ^1^^ ^ _ 3|_^j ^ 

u,+v] = Mu,), j = -l,0,... ,iV, (6.7) 



u 



3 



where A{u) is given by ( |6.5|) . In particular, the first three cases are given by 

p = l: fiiuj) + fiiuj + 7]) + fs{uj) = , 

p = 2: fi{uj)A{u, +7]) + fi{uj + r])A{u, + 27]) + h{uj) = , 

p = 3: fi{uj)A{uj + r])Huj + 2r]) + fi{uj + r])Huj + 2r])A{uj + 3r]) 

+f,{u,) = 0. (6.8) 

The simplest case p = 1 has recently been analyzed in 0.^ 

The Bethe Ansatz equations may be written in a more explicit form by substituting the 



Ansatz ( |6.5|) into ( |6.?1) , and then using the identity ([4.11| ) to simplify the products. In this 



way, we obtain 

• ■ N 2p{N+2) 



X ) fr sinh((p + l){uj - Uk)) sinh((p + l){uj + Uk)) 

I J^^^ smh.{uj — Uk) sinh(Mj — Uk — rj) smh.{uj + Uk + i]) smh.{uj + Uk) 



^The functions gi{u) and 53(14) in Q are related to fi{u) and f^iu) as follows: 

/i(m) — — sinh^ 2u cosh"*^ u gi{u) , fsiu) = ^ sinh^ 2u sinh^^ u cosh^^ u 53 (w) . 
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+ ffu +r])T] sinh((p + l){uj - Uk + rf)) sinh((p + l){uj + Uk + y)) | 

sinh(wj ~ Uk + rj) smh{uj — Uk) smh.{uj + Uk + 21]) smh{uj + Uk + rj) ) 

+ ... = /3K), j = -l,0,... ,iV. (6.9) 

The poles which occur when j = k are harmless, as they are canceled by corresponding zeros. 

We conclude this section with a brief discussion of the diagonal case, k± = 0. In this 
case, the transfer matrix commutes with the matrix Ai, 

[t{u),M] = 0, (6.10) 

where 

N 1 ^ 

M = --S^, S' = ^Y.^l. (6.11) 

i=l 

The two matrices can therefore be simultaneously diagonalized, 

t(M)|A('")) = A('")(m)|A(™)) , 
M\h^"'^) = mlA^™)). (6.12) 

The asymptotic behavior of the eigenvalues is now given by 

AM(„) ^ _l_e«(2iv+2) (e2(7v-™),+5-?+ + e2n.,-c_+c+) ^ _ _ _ ^ oo . (6.13) 

The eigenvalues are therefore given by 

TV 

A('")(m) = p^"*) Yl sinh(M - Uj) sinh(M + r] + uj) , (6.14) 

j=0 

where 

pM = 2ie-'^ cosh (^_ - ^+ + (A^ - 2m)r7) , (6.15) 
and the roots uj satisfy essentially the same Bethe Ansatz equations ( |6.7| ). 



7 Discussion 



We have demonstrated an approach for solving the open XXZ quantum spin chain with 
nondiagonal boundary terms for r] = In particular, we have proposed exact {p + 1)- 
order functional relations ( |5.1|) for the transfer matrix, and Bethe Ansatz equations (|6.7|) , 
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(3.9) for the corresponding eigenvalues. We emphasize that the function fsiu) appearing in 
these equations involves the coefficients Cp^i which we have computed only up to p = 5. (See 
Eqs. ( |5.6| ), ( [4.24| ), ( ^.25| ).) The complete determination of these coefficients must presumably 
await a more explicit construction of nondiagonal K matrices for arbitrary spin. 

A similar approach can certainly also be applied to the elliptic case (i.e., the open XYZ 
chain with nondiagonal boundary terms), since both a fusion hierarchy and a truncation 
identity can also be obtained for this case. Moreover, we expect that a similar approach can 
be applied to spin chains constructed with R and K matrices associated with any affine Lie 
algebra. 

An important feature of our Bethe Ansatz equations is that, as in the case of the closed 
(periodic boundary conditions) XYZ chain, there is a fixed number of roots for all the 
eigenvalues. Because of this fact, it is not easy to compare our results in the diagonal limit 
{k± —>■ 0) with those of |Q, |l|]. A similar difficulty arises ||25[ when comparing the Bethe 
Ansatz equations of the closed XYZ chain in the trigonometric limit with the "conventional" 
Bethe Ansatz equations of the closed XXZ chain. 

It should be interesting to use our results to determine, for low values of p, physical 
properties (thermodynamics, etc.) of the spin chains and of the associated quantum field 
theories. It would also be interesting to try to further simplify our system of Bethe Ansatz 
equations, and perhaps to generalize to the case of generic values of the bulk anisotropy. 
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A Similarity transformations 

We discuss in Section | certain similarity transformations on the fused R and K matrices. 
Here we give explicit expressions for the matrices A , B ,C corresponding to these similarity 
transformations for the first three cases, namely j = 1 , | , 2. 
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For j = 1, 



A 



/ 1 \ 

i ^ 
1 

Vo i -1 0/ 



(A.l) 



B — diag{a ,1 ,a ,1) , C — diag{a ,1 ^1 ,1) , a — {[2\q) ^ 



(A.2) 



For j = |, 



/ 1 














1 

3 











I 
1 

3 









-i i 



1 

3 



1 

2 



1 

3 





1 

3 









1 

3 





- 







1 

3 







1 _2 
3 3 










o\ 





1 







(A.3) 



B — diag{a ,1,1, a, 1,1, 1,1), C — diag{a , 1 , . . . , 1) , a — ([S]^) ^ 



(A.4) 



For j = 2, 

/I 

















vo 




1 

4 






1 

4 
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4 
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6 
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1 

6 
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4 








1 

6 
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6 
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3 








1 

3 





1 
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1 

6 














1 

3 





1 

6 


1 

6 





1 

2 














1 

4 


1 

4 




















1 

4 


1 

4 








1 

4 






1 

4 





1 

3 






1 

2 




1 

4 






1 

4 





1 

3 





_ 1 

2 































1 


1 





1 


6 


6 


6 








1 

4 





























1 


1 





1 


6 


6 


6 








1 

4 

















1 


1 





1 


6 


6 


3 








1 

3 





1 

6 


1 

6 





1 

3 














1 

4 


1 

4 














1 

2 





1 

4 


1 

4 














1 

4 






1 

4 



















1 

4 





_3 
4 












o\ 





1 











0/ 



, (A.5) 



18 



B = diag{a ,l,6,l,a,l,...,l), C = diag{a ,1,6,1,... , 1) 



1 



[2], 



2 



a = ([4],)-^, ^=(j^) • (A-6) 



Note added: 



Boundary quantum group symmetry ^ can be used to determine K matrices 



for arbitrary spin. In particular, the coefficients Cp^i appearing in Eq. (|4.24| ) are given by 



1-2 



k=0 



Thus, for p ^ oo, Cp^i ^. 
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